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ABSTRACT 


Understanding  the  vertical  variations  of  temperature  and  humidity  in  the  marine 
atmospheric  surface  layer  (MASL)  is  extremely  important  for  naval  and  civilian 
applications.  In  particular,  such  variations  affect  the  propagation  of  electromagnetic 
waves  (EM)  by  forming  an  evaporation  duct.  However,  direct  measurements  of  these 
profiles  have  been  difficult  from  a  large  ship  because  of  the  disturbance  introduced  by  the 
platform.  In  this  thesis,  the  design,  deployment,  and  initial  data  analyses  of  a  marine 
atmospheric  profiling  system  (MAPS)  is  introduced.  The  MAPS  is  developed  as  part  of 
the  Coupled  Air  Sea  Process  and  EM  ducting  Research  (CASPER)  project.  It  is  capable 
of  making  repeated  measurements  of  the  lowest  tens  of  meters  of  the  MASE  from  a  small 
Rigid  Hull  Inflatable  Boat  (RHIB),  or  a  small  work  boat,  equipped  with  a  tethered 
profiling  system  and  a  small  meteorological  mast.  Eor  each  profiling  set  at  a  given 
location,  10-15  profiles  were  made  to  allow  sufficient  samples  to  derive  the  mean  profile. 
This  thesis  discusses  the  methods  for  controlling  data  quality  and  obtaining  the  mean 
profiles  from  the  scattered  profiling  data.  Evaporation  duct  height  and  strength  are 
derived  and  compared  to  those  generated  from  an  evaporation  duct  model  using  various 
input  from  measurements. 
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I.  INTRODUCTION 


A.  THESIS  OBJECTIVES 

New  insights  and  future  developments  in  quantifying  boundary  layer  refraction 
profiles  and  therefore  improving  propagation  predictions  will  heavily  depend  on 
unraveling  atmosphere-upper  ocean  processes  on  multiple  scales  using  novel 
measurement  capabilities.  To  understand  the  evaporative  duct,  a  direct  measurement  of 
surface  layer  thermodynamic  profile  with  high  vertical  resolution  is  needed.  This  has  not 
been  done  in  previous  field  efforts.  Unlike  measurements  on  a  tower  or  mast  on  land, 
surface  layer  measurements  over  the  sea  are  difficult  to  obtain  at  multiple  levels.  The 
objective  of  this  thesis  will  be  to  examine  the  feasibility  of  sampling  the  lowest  few  tens 
of  meters  of  marine  atmospheric  surface  layer  (MASL)  in  a  minimally  disturbed 
environment.  As  part  of  the  Coupled  Air-Sea  Processes  and  EM  ducting  Research 
(CASPER)  project,  the  Naval  Postgraduate  School  (NPS)  Meteorology  Department  has 
developed  a  Marine  Atmospheric  Profiling  System  (MAPS)  to  sample  the 
thermodynamic  profiles  in  the  atmospheric  surface  layer.  The  dataset  from  this  system 
will  fill  the  void  of  MASE  profile  measurements  over  the  ocean.  The  system  is  designed 
to  make  profiling  measurements  with  multiple  up/downs  using  an  instrumented  tethered 
balloon  to  increase  the  number  of  samples  at  any  given  altitude  to  provide  high  statistical 
significance.  This  thesis  work  will  demonstrate  the  feasibility  of  the  profiling  system  by 
analyses  of  the  data  quality,  methods  of  deriving  MASE  mean  temperature  and  humidity 
profiles  from  the  raw  measurements,  and  the  derivation  of  evaporative  duct  properties 
based  on  the  measurements.  Eurthermore,  the  dataset  will  be  used  to  evaluate  evaporative 
duct  models  such  as  Navy  Atmospheric  Vertical  Surface  Eayer  Model  (NAVSEaM). 

B.  IMPORTANCE  OF  STUDY 

Correctly  quantifying  the  characteristics  of  the  surface  layer  environment  is 
crucial  to  understanding  low-altitude  Electromagnetic  (EM)  propagation,  particularly  the 
depiction  of  near-surface  moisture  and  temperature  profiles  with  high  vertical  resolution. 
With  the  ultimate  goal  of  improving  evaporative  duct  prediction,  we  use  a  tethered 
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balloon  based  measurement  system  to  obtain  direct  measurements  of  the  surface  layer 
profiles.  This  collection  of  real-time  near  surface  profile  observations  is  essential  for 
evaporative  duct  model  development  and  validation.  EM  propagation  models  also  depend 
on  robust  in  situ  profile  measurements  and/or  highly  accurate  forecasts.  Surface  moisture 
and  temperature  gradients  affect  the  height  of  the  evaporative  duct  that  can  act  as  a  wave 
guide  for  high  frequencies  (Edson  et  al.  1999).  EM  propagation  is  very  sensitive  to  the 
vertical  variation  in  the  atmospheric  boundary  layer  (ABE)  and  would  benefit  from 
improvements  of  model  parameterizations  with  additional  data  for  both  research  and 
operational  purposes.  While  near-surface  data  is  key  to  improving  surface  flux 
parameterization  and  understanding  the  mechanisms  that  couple  the  ocean  and 
atmosphere,  it  remains  one  of  the  greatest  challenges  to  obtain  the  data  observationally. 
Ultimately,  accurate  prediction  of  the  surface  layer  profiles  and  fluxes  is  the  key  to 
weather  prediction  and  EM  prediction. 

Bulk  aerodynamic  surface  flux  parameterization  schemes  are  needed  to  explain 
the  exchanges  of  mass  and  energy  across  the  air-sea  interface.  These  are  largely  built 
upon  the  Monin-Obukhov  similarity  theory  (MOST)  developed  in  1954  (Eiu  et  al.  1979). 
Eairall  et  al.  (1996b;  2003)  further  modified  the  bulk  surface  flux  system  by  using  data 
obtained  during  the  Tropical  Ocean  Global  Atmosphere  (TOGA)  Coupled  Ocean- 
Atmosphere  Response  Experiment  (COARE)  that  resulted  in  what  is  known  as  the 
COARE  bulk  flux  algorithm.  COARE  is  now  the  most  widely  used  flux  algorithm  in 
most  of  the  mesoscale  and  global  scale  forecast  models.  Eimitations  to  the  COARE 
algorithm  do  exist,  however,  partly  due  to  the  inherent  assumptions  of  MOST  (e.g., 
Andreas  et  al.  2014).  Thus,  the  need  to  collect  high-quality  data  at  the  air-sea  interface 
continues  in  order  to  further  improve  surface  flux  parameterization  algorithms  and  hence 
forecast  models. 

C.  NAVAL  APPLICATION 

Accurately  characterizing  the  spatial  and  temporal  variability  of  refractivity  is 
crucial  to  many  important  Navy  and  civilian  applications.  The  ultimate  goal  of  the  project 
is  to  enhance  the  Navy’s  capability  of  predicting  the  performance  of  radar  and 
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communications  systems  for  tactical  applications.  With  the  removal  of  the  rawinsonde 
program  from  the  U.S.  Navy  fleet  in  2011,  the  accuracy  of  EM  propagation  prediction  for 
surface  and  airborne  radar  and  communication  operations  depends  primarily  on  the 
modeled  environmental  forecasts  that  feed  the  propagation  models.  “EM  propagation  is 
sensitive  to  even  slight  changes  in  the  ABE  temperature  and  moisture  gradients”  (Babin 
1997).  Near  surface  measurements  are  critical  to  determining  and  accurately  forecasting 
these  gradients  and  the  structure  of  the  ABE.  Advanced  Refractive  Effects  Prediction 
System  (AREPS)  and  Target  Acquisition  Weapons  Software  (TAWS),  in  particular, 
require  accurate  ABE  profile  inputs  to  effectively  predict  EM  propagation,  radar  ranges, 
and  weapon  sensor  effectiveness. 

The  Navy  is  committed  to  achieving  superiority  of  the  electromagnetic  spectrum 
(EMS).  The  2015  Electromagnetic  Maneuver  Warfare  (EMW)  Campaign  Plan  “is  the 
Navy’s  warfighting  approach  to  gain  decisive  military  advantage  in  the  EMS  to  enable 
freedom  of  action  across  all  Navy  mission  areas.”  Commander  Navy  Meteorology  and 
Oceanography  Command,  Rear  Admiral  Gallaudet  has  outlined  a  strategy  to  advance 
Navy’s  electromagnetic  warfare  capabilities  in  a  recent  article.  One  of  the  goals  of  this 
strategy  is  to  improve  Naval  Oceanography’s  Environmental  and  Prediction  Capabilities 
by  advancing  our  environmental  sensing  capabilities: 

The  need  for  additional  sensing  is  most  critical  in  the  lowest  portion  of  the 
atmosphere  where  the  evaporative  duct  forms  and  impacts  the  propagation 
of  signals.  We  will  develop,  evaluate  and  transition  autonomous/ 
unmanned  platforms  and  sensors  that  enable  persistent,  physical 
battlespace  awareness,  potentially  through  a  new  Eittoral  Battlespace 
increment.  We  will  improve  the  vertical  resolution  of  our  observations  in 
the  lowest  portion  of  the  atmosphere  where  impacts  on  EMW  propagation 
is  the  highest.  (Gallaudet  2016) 
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II.  BACKGROUND 


A.  ATMOSPHERIC  REFRACTION 

“Atmospheric  refractive  conditions  can  significantly  affect  the  performance  of 
shipboard  radars  and  communications  at  sea  and  near  shore”  (Battan  1973). 
Characterizing  the  spatial  and  temporal  variability  of  refractivity  is  thus  crucial  to  many 
important  Navy  and  civilian  applications.  When  electromagnetic  radiation  travels  through 
the  atmosphere,  they  follow  curved  instead  of  straight  paths  such  as  in  the  outer  space 
without  the  atmosphere.  The  bending  of  radio  waves  is  usually  quantified  by  the  index  of 
refraction  defined  as  the  ratio  between  the  speed  of  light  in  a  vacuum  over  that  in  a 
medium.  The  index  of  refraction  is  affected  by  temperature,  pressure  and  humidity,  the 
latter  being  the  most  important,  as  described  in  Equation  (1)  (Bean  and  Dutton  1968): 

A  =  (n-1)  X  10^  =  — (^p+— j  (1) 

where  N  is  refractivity,  n  is  the  index  of  refraction,  T  is  the  atmospheric  temperature  (K), 
p  is  the  total  atmospheric  pressure  (hPa),  and  e  is  the  water  vapor  pressure  (hPa).  To 
consider  propagation  over  the  earth  with  a  curvature,  a  modified  refraction,  M,  is  defined 
so  that  one  can  considered  the  earth  as  a  hypothetically  flat  surface: 

M  =  A +  ^^(n- 1)  X  10^  =  iV  +  0.157Z  (2) 

where  Vg  is  the  earth’s  radius  in  meters  and  z  is  altitude  in  meters  (Bean  and  Dutton 
1968). 


The  refractivity  profile  in  the  atmosphere  determines  the  curvature  of  the  ray 
depends  on  the  rate  of  change  of  the  refractive  index  with  height. 

Figure  1  illustrates  the  refractivity  propagation  categories  and  the  ray-paths  of 
horizontally  transmitted  rays.  These  categories  are  based  on  the  slopes  of  N  (dN/dz)  or 
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M  {dM/dz)  for  sub-refraction,  standard  refraction,  super-refraction,  and  trapping  (Turton 
1988).  In  a  paper  written  by  Turton  (1988),  he  explains  that  when  the  refractivity  gradient 
is  greater  than  157  M  units  km'\  the  EM  wave  will  bend  upwards.  These  conditions  are 
said  to  be  sub -refractive  and  result  in  short  detection  ranges.  Standard  refractive 
conditions  is  when  the  modified  refractivity  gradient  is  between  78  and  157  M  units  km'\ 
These  are  considered  normal  propagation  in  contrast  to  the  rest  of  propagation  categories, 
which  are  normally  referred  to  as  abnormal  propagation  (AP).  When  the  M  gradient  is 
less  than  78  M  units  km'\  the  conditions  are  defined  as  super-refractive  where  the  EM 
rays  bend  downward  compared  to  standard  refraction.  If  the  M  gradient  is  less  than  zero, 
the  ray’s  curvature  exceeds  that  of  the  earth  it  becomes  trapped  in  a  quasi-horizontal  layer 
called  a  duct.  In  super  refractive  conditions,  and  particularly  in  ducting  conditions,  radar 
detection  ranges  become  significantly  increased  if  the  source  is  within  the  duct  (Turton 
1988). 


<JN/d/  =  0 
dM/dz  =  157 


Sub-refraclion  ^ '' 


Radar 


Eigure  1 .  Propagation  categories  and  the  ray-paths  of  horizontally  transmitted 
rays.  Categories  are  based  on  the  slopes  of  N  (dN/dz)  and  M  (dM/dz)  for 
sub-refraction,  standard  refraction,  super-refraction,  and  trapping. 
Source:  Turton  (1988). 
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B.  EVAPORATIVE  DUCT 


Evaporation  ducts  are  ever  present  over  the  ocean,  and  are  associated  with  strong 
vertical  water  vapor  gradients  near  the  sea  surface  (Babin  et  al.  1997).  This  type  of  duct 
occurs  in  the  MASL  as  a  result  of  sea  surface  evaporation  and  this  evaporation  is  enough 
to  cause  the  bending  of  electro-magnetic  radiation  so  that  it  becomes  trapped  in  this  area 
(i.e.,  ducting)  (Babin  et  al.  1997).  The  evaporation  duct  is  hence  a  special  type  of  surface 
duct  with  the  trapping  layer  extending  to  the  surface.  In  identifying  the  trapping  layers  of 
the  troposphere,  the  modified  refractivity  expressed  in  Equation  (2)  is  most  often  used  to 
characterize  the  evaporation  duct  properties.  Traditionally,  according  to  a  paper  written 
by  Babin  and  Dockery  (2001),  “the  evaporation  duct  has  been  characterized  by 
determining  only  the  height  of  the  duct  as  defined  by  the  altitude  at  which  dM/dz  =  0. 
However,  it  is  the  slope  dM/dz  that  is  used  in  propagation  models,  and  a  given  duct 
height  may  result  from  a  variety  of  different  shapes  of  M  profiles.  As  Equation  (1)  and 
(2)  indicate,  these  M-profile  curvatures  are  primarily  affected  by  vertical  temperature  and 
moisture  profiles  and  hence  by  atmospheric  static  stability.”  Gehman  points  out  in  a 
paper  written  in  2000,  “the  actual  slope  is  more  important  in  uniquely  determining 
electromagnetic  propagation  effects  than  the  duct  height  alone.”  Duct  strength  can  be 
identified  by  examination  of  the  vertical  profile  of  the  modified  refractivity.  Eigure  2 
represents  an  M  profile  under  evaporative  duct  conditions  and  illustrates  how  the  duct 
depth  and  strength  are  determined.  The  duct  layer  is  from  the  altitude  of  M’s  local 
minimum  to  the  surface,  which  is  also  and  the  trapping  layer. 
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Modified  refractivity 

Figure  2.  Definition  of  evaporative  duct  height  and  duct  strength. 

Source:  (Turton  et  al.)  1988. 

Electromagnetic  propagation  is  usually  affected  by  evaporative  ducts  within  the 
first  40  m  of  the  oceanic  surface.  Babin  and  Dockery  (2001)  explain  that  “they  are  of 
particular  importance  to  naval  operations:  ducting  can  significantly  enhance  the  range  of 
radio  communications  and  the  radar  detection  of  targets.”  Ducting  can  dramatically 
improve  the  distance  of  radio  communications  and  radar  detection  of  targets  (Babin  and 
Dockery  2001).  On  the  other  hand,  reduced  radio  communications  and  radar  coverage 
(radar  holes)  can  occur  just  outside  the  trapping  layer  due  to  the  focus  of  energy  inside 
the  duct.  For  this  reason,  ducting  is  of  significance  importance  to  naval  operations  (Babin 
and  Dockery  2001).  Figure  3  depicts  the  radio  propagation  path  in  case  of  an  evaporative 
duct. 


8 


Figure  3.  Radio  propagation  path  in  case  of  evaporative  duct. 
Source:  Turton  et  al.  (1988). 


C.  MARINE  ATMOSPHERIC  SURFACE  LAYER  FLUX-  PROFILE 

RELATIONSHIP  AND  MONIN-OBUKHOV  SIMILARITY  THEORY 

In  a  comprehensive  paper  written  by  Edson  et  al.  (2003),  it  states  that  at  the  air 
ocean  interface,  a  significant  amount  of  energy,  momentum,  and  mass  is  transferred 
between  the  atmosphere  and  the  ocean  in  the  lower  boundary  layer.  These  exchanges  are 
carried  out  in  the  form  of  water  vapor,  momentum,  and  heat.  These  variables  are 
measured  by  momentum,  latent  and  sensible  heat  fluxes.  The  transfer  rate  of  these  fluxes 
can  have  a  significant  effect  on  the  structure  of  the  MASL  and  the  ocean  mixed  layer 
(Edson  et  al.  2003).  These  scholars  contend  that  one  of  the  most  challenging  areas  of 
environmental  forecasting  is  representing  the  interaction  of  the  air-sea  processes  and  their 
effects  in  environmental  diagnostic  and  prognostic  models.  In  the  surface  layer,  turbulent 
eddies  transfer  momentum,  latent  and  sensible  heat  fluxes  throughout  the  surface  layer. 
Evaporation  is  the  main  process  that  adds  water  vapor  to  the  atmosphere  in  the  form  of 
water  vapor.  Evaporation  rate  is  influenced  by  the  following  factors:  air  temperature, 
wind  speed,  humidity,  sea  surface  temperature,  and  sea  state  (Edson  et  al.  2003). 

The  lowest  10%  of  the  boundary  layer  is  commonly  thought  of  as  the  surface 
layer.  In  this  layer,  flux  values  typically  do  not  differ  by  more  than  10%  (Stull  1988).  Eor 
this  reason,  theory’s  based  on  physics  in  the  surface  layer  can  be  utilized  to  calculate 
fluxes  for  the  entire  surface  layer  from  the  fluxes  at  one  level. 

In  the  atmospheric  surface  layer,  the  primary  sources  for  generating  turbulence 
are  buoyancy  and  wind  shear.  Because  the  evaporation  duct  occurs  within  the  MASE, 
theories  on  the  physics  of  the  surface  layer,  such  as  Monin-Obukhov  Similarity  Theory 
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(MOST)  are  used  for  evaporative  duct  modeling.  MOST  is  based  on  scaling  analysis  to 
derive  the  relationship  between  the  mean  wind,  temperature,  and  specific  humidity 
profiles  to  surface  layer  turbulent  fluxes  (Edson  et  al.  2003).  The  resultant  non- 
dimensional  gradients  can  be  expressed  as: 


KZ  du  _  fz\ 

Ut  dz  VL/ 


KzdO  f  z\ 

e.  dz  ~  \l) 


Kzdq  f  z\ 

dz  ~  \l) 


(3) 

(4) 

(5) 


where  von  Karman’s  constant  k  ~  0.35  —  0.40  derived  from  measurements,  L  is  the 
Monin-Obukhov  length:  L  =  — — —  ,  g  is  gravity,  z  is  height,  and  f-)  and  f-)  are 

KQWfOf  \L/  \L/ 

non-dimensional  U  functions  that  account  for  the  effects  of  thermal  stability  of  the 
surface  layer.  In  addition,  u*  is  the  frictional  velocity,  q*  is  specific  humidity  scale  and  0* 
is  the  temperature  scale.  According  to  Edson,  these  are  considered  the  scaling  parameters 
of  the  surface  layers.  The  left  hand  side  of  the  equations  represents  the  non-dimensional 
vertical  gradient  of  mean  wind,  mean  potential  temperature,  and  mean  specific  humidity. 
The  right  hand  side  of  the  equations  shows  that  these  mean  gradients  are  functions  of 
height  non-dimensionalized  by  L,  signifying  the  effects  of  thermal  stability.  These 
universal  functions  have  been  empirically  derived  from  measurements  in  the  field. 
Businger  et  al.  (1971)  presented  their  commonly  used  land-based  universal  functions  in 
the  famous  1968  Kansas  experiment: 


10 


unstable 


{1-15-)-^'" 

L 


(1  +  4.7-) 
L 


z  _Q  neutral 

L~ 

z  „  stable 

—  >0 
L 


0.74(1-9-) 

L 


^  ^  -1/2 


0.74 


(0.74-1-4.7^) 


z  ^  unstable 
—  <0 
L 


z  _Q  neutral 
L~ 


z  „  stable 
—  >0 
L 


The  MOST  non-dimensional  functions  are  responsible  for  the  connection  between 
mean  profiles  and  turbulence  measurements.  This  connection  developed  by  Monin  and 
Obukhov  provides  two  main  functions.  First,  MOST  provides  a  method  to  calculate 
surface  layer  turbulent  fluxes  using  mean  profile  measurements.  This  process  is 
particularly  important  when  direct  sampling  of  turbulence  is  not  possible.  MOST  also  lets 
researchers  calculate  mean  wind  at  a  level  at  which  a  direct  observations  was  not  taken. 
For  example,  many  flux  parameterization  equations  need  mean  winds  values  at  10  m.  By 
using  the  non-dimensional  gradient,  known  mean  values  at  one  level,  and  surface  fluxes, 
researchers  can  compute  the  wind  at  10  m.  If  observations  are  obtained  at  only  one  level 
and  flux  is  not  known,  both  the  gradient  and  the  surface  turbulent  fluxes  can  be  derived 
by  using  the  surface  roughness  height  and  assumed  mean  properties  at  the  roughness 
length  (zo),  defined  as  the  height  at  which  mean  wind  is  zero  (Chamock  1955).  Also,  the 
air  temperature  at  the  surface  roughness  height  is  set  equal  to  SST,  and  the  relative 

humidity  over  salt  water  at  the  surface  roughness  height  is  taken  as  98%. 
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The  stability  function  in  MOST  was  developed  from  measurements  over  land, 
which  brings  into  question  its  validity  over  the  ocean.  The  fluid  interface  over  the  ocean 
whose  roughness  changes  with  the  wind  generated  wave  and  swells  from  distant  storms 
makes  the  MASL  much  more  complicated  than  its  counterpart  over  land.  Additionally, 
MOST  assumes  horizontally  homogenous  conditions.  It  was  found  (Tellado  2013)  that 
significant  inhomogeneity  always  exists  over  the  coastal  ocean  near  Monterey  Bay.  The 
homogenous  assumption  thus  may  not  hold  true  near  the  coast  or  close  to  a  frontal 
boundary. 

D.  EVAPORATIVE  DUCT  MODELS 

According  to  paper  written  by  Babin  (1997),  “The  general  approach  for  all 
evaporation  duct  models  (EDM)  involves  finding  an  expression  for  the  vertical 
refractivity  gradient  in  terms  of  atmospheric  variables.”  Since  the  evaporative  duct 
occurs  in  the  surface  layer,  MOST  is  typically  used  to  derive  the  mean  vertical  profile. 
EDM’s  calculate  the  modified  refractivity  gradient,  from  this,  evaporative  duct  properties 
such  as  duct  strength  and  height  can  be  obtained  by  evaluating  the  profile  (Cherrett 
2015).  A  comprehensive  summary  of  this  approach  is  provided  in  Babin  (1997).  There 
are  currently  two  types  of  evaporation  duct  models:  the  potential  refractivity  model  and 
the  EKB-based  models. 

Erom  1978  until  2012,  the  Paulus-Jeske  (PJ)  evaporation  duct  model  (Jeske  1973; 
Paulus  1984,  1985,  1989)  was  the  U.S.  Navy’s  most  widely  used  evaporation  duct  model 
(Babin  1997).  It  was  incorporated  into  Integrated  Refractive  Effects  Prediction  System 
(IREPS)  and  then  AREPS  Tactical  Decision  Aids  (TDA).  The  PJ  model  uses  the  potential 
refractivity  quantity,  which  calculates  refractivity  index  using  “potential  temperature, 
potential  water  vapor  pressure,  and  1000  hPa  atmospheric  pressure  as  opposed  to  air 
temperature,  water  vapor  pressure,  and  surface  pressure”  (Jeske  1973).  PJ  model  allocates 
values  of  wind  speed,  relative  humidity,  and  air  temperature  to  a  height  of  6  m,  regardless 
“of  the  actual  observation  height.  SST  is  also  used  and  surface  pressure  is  assigned  a 
constant  value  of  1000  hPa”  (Paulus  1984).  An  additional  parameter  distinctive  to  the  PJ 
model  is  that  instead  of  using  the  typical  -0.157  critical  gradient  for  potential  refractivity 
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it  uses  a  critical  potential  refractivity  gradient  of  -0.125  to  determine  ducting  (Babin 
1997).  Cook  and  Burk  (1992)  indicated  certain  inadequacy  of  the  PJ  model.  They 
revealed  that  “when  properly  non-dimensionalized  the  vertical  gradient  of  potential 
refractivity  was  not  a  single  universal  function  of  z/L  and  that  potential  refractivity  in 
stable  conditions  did  not  obey  MOST”  (Babin  1997).  This  implied  that  for  stable 
conditions,  the  principle  theory  for  this  model  in  using  potential  refractivity  was 
inappropriate. 

Liu  et  al.  (1979)  developed  a  marine  atmospheric  surface  layer  model  that  is 
referred  to  as  Liu,  Katsaros,  and  Businger  (LKB)  and  used  it  for  deriving  air-sea 
exchanges  of  heat,  moisture,  and  momentum.  The  LKB  model  included  the  interfacial 
molecular  effects  at  the  sea  surface  and  matched  the  mean  wind  and  scalar  profiles  from 
the  MOST  theory  with  those  in  the  molecular  sublayer.  Their  approach  was  used  to 
represent  surface  layer  turbulence  fluxes,  which  leads  to  the  MOST-based  surface  flux 
parameterizations.  The  Tropical  Ocean  and  Global  Atmosphere  Coupled  Ocean- 
Atmosphere  Response  Experiment  (TOGA  COARE)  (Eairall  et  al.  1996)  project  later 
refined  the  EKB  MASE  models  by  incorporating  the  results  of  more  than  10,000  hours  of 
atmospheric  and  oceanic  measurements  from  buoys,  ships  and  aircraft.  Their  efforts  led 
to  the  COARE  surface  flux  algorithm  that  is  currently  the  most  widely  used  surface  flux 
parameterization  scheme  extensively  used  in  various  mesoscale  and  global  scale  models 
(Eairall  et  al.  1996,  2003). 

There  are  several  EKB-based  evaporation  duct  models;  the  most  significant 
difference  among  them  is  the  use  of  different  stability  function,!/;  (z/L)  to  “account  for 
deviations  from  neutral  stability  in  the  atmospheric  surface  layer  and  are  integrated  forms 
of  the  dimensionless  gradient  functions”  (Blackadar  1993).  The  most  notable  EKB-based 
model  is  the  Naval  Postgraduate  School  (NPS)  evaporation  duct  model  described  in 
Erederickson  et  al.  (2000).  In  2012,  The  NPS  model  was  developed  and  integrated  into 
AREPS  and  is  now  called  the  NAVSEaM  (Navy  Atmospheric  Vertical  Surface  Eayer 
Model). 
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E.  COUPLED  AIR-SEA  PROCESSES  AND  ELECTROMAGNETIC 

DUCTING  RESEARCH 

The  Coupled  Air-Sea  Processes  and  EM  ducting  Research  (CASPER),  is 
sponsored  by  the  Office  of  Navy  Research  (ONR)  2014  Multi-disciplinary  University 
Research  Initiative  (MURI)  to  address  overarching  knowledge  gaps  related  to 
electromagnetic  wave  propagation  in  coastal  MABE.  The  objective  of  CASPER  is  to 
fully  characterize  the  MABE  as  an  EM  propagation  environment.  There  are  three 
components  to  the  CASPER  project:  theoretical  developments,  field  program,  and 
numerical  modeling  efforts.  This  thesis  is  centered  on  the  field  program  that  focuses  on 
employing  new  environmental  measurement  techniques  and  novel  sampling  strategies  to 
obtain  a  comprehensive  and  cohesive  dataset  to  address  air-sea  interaction  processes  that 
affect  EM  propagation  and  for  extensive  model  evaluation  and  testing  and  new 
measurement  capabilities.  "The  field  components  have  two  main  campaigns:  CASPER- 
East  (Duck,  NC)  and  CASPER-West  (Southern  California).  With  these  two  operations, 
CASPER  plans  to  fully  characterize  the  marine  atmospheric  boundary  layer  (MABE)  as 
an  electromagnetic  (EM)  propagation  environment.  The  emphasis  will  be  on  spatial  and 
temporal  heterogeneities  and  surface  wave/swell  effects”  (Wang  et  ah). 

Prior  to  onset  of  Casper  East,  the  CASPER  pilot  experiment  was  conducted 
offshore  and  at  the  shoreline  of  Moss  Eanding,  CA.  to  test  out  a  few  of  the  key  platforms 
and  sensors.  CASPER  pilot  took  place  from  April  20,  2015  to  May  2,  2015,  with 
CASPER-East  being  conducted  off  the  coast  of  Duck,  North  Carolina  from  October  10, 
to  November  6,  2015.  All  data  used  in  this  study  were  collected  during  CASPER  Pilot 
and  CASPER  East. 
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III.  NPS  MASL  PROFILING  SYSTEM 


A.  SYSTEM  AND  SENSOR  INFORMATION 

We  have  developed  a  profiling  system  to  sample  the  MASL  in  an  undisturbed 
environment.  The  NPS  Marine  Atmospheric  Surface  Layer  Profiling  system  (NPS 
MAPS)  includes  two  components:  the  profiling  and  ship  mast  components.  The  profiling 
component  intends  to  sample  the  variation  of  temperature  and  relative  humidity  with 
height  (T/RH  profiles)  with  multiple  ascends  and  descends  of  the  probes,  while  the  ship 
mast  component  provide  auxiliary  data  collected  continuously  at  a  fixed  level.  The 
auxiliary  data  includes  wind  speed  and  direction,  pressure,  GPS  position,  and  sea  surface 
temperature  (SST).  Temperature  and  dew  point  temperature  are  also  measured  from  the 
mast.  The  data  from  both  components  are  used  to  derive  refractivity  profiles  used  to 
identify  the  characteristics  of  evaporation  duct  and  to  make  diagnostic  model  calculations 
of  the  evaporation  duct. 

A  prototype  of  the  profiling  system  was  tested  in  CASPER  pilot  experiment  in 
April/May  of  2015.  An  improved  system  was  used  in  CASPER  East  field  campaign.  In 
the  following  description,  both  the  prototype  system  and  the  improved  system  will  be 
described. 

I.  The  Profiling  Component 

The  profiling  component  is  mainly  composed  of  a  radiosonde,  a  tethered  balloon, 
and  a  radiosonde  receiver/display  system.  It  also  includes  a  reel  that  controls  the  up  and 
down  movement  of  the  tethered  balloon.  The  radiosonde  is  the  iMet-l-ABX  sonde  that 
measures  temperature,  relative  humidity,  pressure,  GPS  time  and  locations.  Since  a 
normal  rawinsonde  on  a  free  flying  balloon  obtains  wind  from  the  time  and  position  of 
the  sonde,  the  tethered  sonde  cannot  provide  wind  measurements.  In  Table  1  below, 
different  specifications  of  the  profiling  component  are  compared  between  the  CASPER 
pilot  and  CASPER  East  deployments.  CASPER  East  represents  significant  improvements 
from  the  prototype  in  CASPER  pilot.  The  real-time  transmission  and  display  of  the 
profile  is  a  major  improvement  because  problematic  probes  or  some  sampling  issues  can 
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be  identified  and  corrected  on  the  spot  to  ensure  measurement  success.  The  use  of  the 
electric  reel  made  it  possible  for  multiple  ascends/descends  without  extreme  fatigue  of 
the  operator,  especially  in  windy  conditions.  Using  a  smaller  sized  balloon  and  thinner 
tetherline  also  helped  to  reduce  the  excess  pulling  of  the  balloon. 


Table  1.  Comparison  of  the  profiling  components  used  in  CASPER  Pilot  and 

CASPER  East. 


Specifications 

CASPER  Pilot 

CASPER  East 

Sonde 

iMet-l-ABX 
modified  to  record 
data  on  microchips 

iMet-l-ABX 

Data  collection 

Self-recording 

Transmitted  and 
received  in  real¬ 
time  using  radio 
receiver 

Data  display 

None 

Real-time  display 
on  tablet 

Sonde  attachment 

Separately  tethered 
free  dangling  -5  m 
below  the  balloon 

Attached  directly 
to  the  tetherline  of 
the  balloon  ~3  m 
below  the  balloon 

Balloon 

3.3  m^  Helikite 

2.0  m^  Helikite 

Balloon  motion  control 

manual 

Electric  fishing  reel 

Power  need 

None 

Marine  battery  to 
power  fishing  reel 

In  both  CASPER  pilot  and  East,  we  used  the  Helikite  as  the  tethered  lifting 
system.  The  Helikite  aerostats  are  kite-balloon  hybrid  aerostat.  In  low  wind  conditions, 
the  lift  is  mainly  from  the  lighter-than-air  balloon,  while  in  wind  conditions,  the  kite 
produces  significant  lift.  The  major  advantage  of  the  Helikite  is  its  stability  compared  to 
normal  tethered  balloons  and  kites,  especially  near  the  surface  where  there  is  significant 
turbulence.  In  conditions  with  some  wind,  the  Helikites  fly  at  about  45°  angle.  The  2  m^ 
Helikite  used  in  CASPER  East  is  9  ft  x  7  ft  in  size  (length  x  width)  and  provides  0.8  kg 
of  lift  in  no  wind  conditions  and  ~4  kg  of  lift  in  wind  speed  of  ~6.7  ms’^  (15  mph).  When 
fully  deployed,  the  radiosonde  flies  about  50  m  above  the  waterline  and  when  fully  reeled 

in  the  radiosonde  lies  below  1  m  from  the  ocean  surface. 
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2. 


Sensors  on  Mast 


The  small  mast  on  the  ship  hosts  a  suite  of  sensors  that  provide  auxiliary  data  for 
ED  property  retrieval  and/or  EM  model  input.  These  sensors  include  a  gyro  stabilized 
electronic  compass  which  provides  accurate  heading,  pitch,  and  roll  in  various  dynamic 
conditions,  a  Vaisala  Weather  Transmitter  WXT520  that  measures  barometric  pressure, 
humidity,  temperature,  wind  speed  and  direction.  Sea  surface  temperature  was  measured 
by  a  thermistor  located  about  0.02  m  below  the  sea  surface.  The  sensors  and  the  data 
acquisition  system  were  all  powered  by  a  12V  marine  battery. 

In  Casper  Pilot,  A  20  ft  RHIB  with  a  bow  walk  that  protrudes  out  over  the  water. 
The  met  mast  was  setup  on  this  bow  walk  away  from  operator  activities  to  avoid  flow 
contamination  (Eigure  4).  The  mast  sensors  were  at  a  height  of  about  2  m  for  all  data 
collection  during  the  pilot  experiment.  The  setup  in  CASPER  East  is  shown  in  Eigure  5, 
where  the  work  boat  of  the  R/V  Hugh  Sharp  was  used.  The  mast  was  attached  to  the 
metal  structure  on  the  opposite  side  of  the  tethered  balloon  attachments  and  operation. 

The  small  boat  set-up  was  significantly  improved  for  CASPER-East  deployment. 
An  electric  fishing  reel  replaced  the  manual  reel  used  in  the  pilot  experiment.  This 
improved  the  operation  of  the  tethered  balloon  to  allow  for  multiple  vertical  profiles 
during  each  profiling  session  without  excessive  operator  fatigue.  The  MET  mast  was  also 
moved  further  away  from  the  operators  to  avoid  flow  contamination.  Easily,  a  smaller 
balloon  (2-m^)  replaced  the  larger  Helikite  (3.3-m^)  used  in  CASPER- Pilot,  which 
eliminated  the  excessive  pull  from  the  tethered  system  and  hence  less  requirement  on  the 
reeling  system. 
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Helikite 


i  Self-recording  sonde 


Figure  4.  Prototype  RHIB-based  tethered  balloon  MAPS  used  in  CASPER 
Pilot.  The  main  components  of  the  system  are  annotated. 
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Figure  5.  NFS  MAPS  deployed  in  CASPER  East.  The  electric  reel  and  the 
sensor  mast  are  enlarged  as  inserts. 

3.  Sensor  Accuracy 

All  sensors  of  the  MAPS,  their  sampled  variables  and  the  corresponding  accuracy 
are  listed  in  Table  2.  Table  2  shows  that  temperature  and  relative  humidity  on  the  mast 
are  sampled  redundantly  from  two  probes.  The  Rotronic  probe  was  added  because  of  its 
better  accuracy  for  T/RH.  Since  the  sonde  was  tethered,  no  wind  information  could  be 
derived  from  the  sonde  measurements.  The  Vaisala  WXT520  provides  mean  wind  from 
measurements  at  the  mast  level.  The  GPS  and  compass  provide  the  platform  locations 
and  can  be  used  to  derive  platform  speed  when  needed.  Due  to  the  need  of  fast  response 
time  for  profiling  measurements,  the  response  time  and  sensor  resolution  of  the  iMet 
sondes  are  also  given  in  Table  2.  Note  that  although  the  sampling  rate  of  the  sonde  is  at  1 
Hz  as  specified  by  the  manufacturer,  given  that  the  normal  ascend/descend  speed  were 
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controlled  at  about  0.5  ms'^  and  the  response  time  of  the  sonde  probes  in  Table  2,  the 
MAPS’  measurements  of  temperature  and  humidity  have  an  effective  vertical  resolution 
of  nominally  1  m. 


Table  2.  Accuracy  of  all  sensors  of  the  NPS  MAPS 


Sensor 

Measured  variables  and  accuracy 

Vaisala  WXT520  Weather  station 

Temperature/Humidity  (±0.3°C/±3%) 

Wind  speed  and  direction  (±0.3  ms'^  and  ±3°) 
Barometric  pressure  (±0.5  hPa) 

Precipitation  (5%) 

Rotronic  HC2-S3 

Temperature  (±0.1°C) 

Relative  humidity  (±0.8%) 

Garmin  GPS16x-HVS 

Position  (<15  m),  GPS  time 

Course  over  ground  and  speed  over  ground 
Magnetic  declination 

TNT  Compass 

Heading  (±3°) 

Pitch/roll  (<1°) 

iMET  rawinsonde 

Pressure  (response  time  <  1.0  sec,  accuracy  0.5  hPa,  resolution 
0.01  hPa) 

Temperature  (response  time  2  sec,  accuracy  0.2°C,  resolution 

0.01  °C) 

Relative  humidity  (response  time  2  sec,  accuracy  5%,  resolution 

<  0.01%) 

Top  layer  water  temperature 

Temperature  (±0.1°C) 

B.  FIELD  DEPLOYMENT 

As  previously  stated,  all  data  used  in  this  study  was  collected  during  CASPER 
Pilot  (Figure  6)  and  CASPER  East  (Figure  7).  These  two  field  campaigns  captured 
ducting  conditions  off  shore  of  California  coast  as  well  as  on  the  East  coast.  The  NPS 
MAPS  cruise  track  in  CASPER  Pilot  was  orientated  in  east- west  direction  and  stretched 
approximately  11  nautical  miles  offshore.  Three  locations  along  the  linear  track  were 
designated  to  obtain  profile  measurements.  In  CASPER  East  MAPS  was  mostly 
orientated  in  west-east  direction  and  measured  approximately  25  nautical  miles  offshore. 
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Latitude  ('N) 


Longitude  (°W) 

Figure  6.  MAPS  deployment  locations  (red  dots)  and  the  ship  (RHIB)  track 
(blue  lines)  during  CASPER  Pilot  offshore  of  Moss  Landing,  California. 
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35.5N 


76W  75.5W  75W  74.5W  74W 


Figure  7.  MAPS  deployment  locations  (red  dots)  and  the  ship  (R/V  Hugh 
Sharp)  track  (blue  lines)  during  CASPER  East  offshore  of  Duck,  North 

Carolina. 
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At  each  profiling  location,  a  minimum  of  seven  up  and  down  vertical  profiles 
(yielding  a  set  of  14  profiles  each  time)  were  obtained  by  releasing  the  balloon  to  a  height 
of  about  50  meters  and  then  reeling  the  sensor  package  down  to  close  to  the  surface.  In 
both  field  campaigns,  the  RHIB/work  boat  was  also  equipped  with  an  X-band  beacon 
transmitter  so  that  RF  propagation  between  the  RHIB/work  boat  and  the  main  ship  can  be 
measured.  Each  set  of  measurements  was  completed  within  approximately  30  minutes 


Table  3 .  CASPER  Pilot  NPS  MAPS  profiling  data 


Deployment  # 

Date 

Start  time 

(UTC) 

End  time 

(UTC) 

Duration 

(Minutes) 

#  of 

up/downs 

1 

20-Apr 

21:57 

22:30 

33 

7 

2 

21-Apr 

18:53 

19:30 

27 

6 

3 

26-Apr 

17:52 

18:22 

30 

6 

4 

26-Apr 

19:02 

19:33 

31 

7 

5 

27-Apr 

17:16 

17:45 

29 

7 

6 

27-Apr 

18:00 

18:32 

32 

7 

7 

27-Apr 

18:30 

18:59 

29 

7 

8 

01-May 

16:19 

16:43 

24 

7 

9 

01-May 

16:53 

17:20 

27 

7 

10 

01-May 

17:24 

17:46 

22 

6 

11 

01-May 

20:20 

20:50 

30 

7 

12 

01-May 

20:53 

21:21 

28 

7 

11 

02-May 

15:52 

16:22 

30 

7 

12 

02-May 

16:28 

17:30 

32 

7 

13 

02-May 

17:45 

18:10 

25 

6 

14 

02-May 

18:34 

19:02 

28 

7 

12 

02-May 

21:06 

21:36 

30 

7 

13 

02-May 

21:30 

21:58 

28 

7 

14 

02-May 

22:00 

22:26 

26 

7 

23 


Table  4.  PER-East  NPS  MAPS  pronling  data 


Deployment#  Date 

1  13-Oct 

2  _ 13-Oct 

3  _ 14-Oct 

4  _ 14-Oct 

5  _ 14-Oct 

6  _ 15-Oct 

7  _ 15-Oct 

8  _ 15-Oct 

9  _ 16-Oct 

W _ 16-Oct 

n _ 16-Oct 

U _ 17-Oct 

U _ 17-Oct 

M _ 17-Oct 

15 _ 17-Oct 

15 _ 20-Oct 

17 _ 20-Oct 

15 _ 20-Oct 

15 _ 21-Oct 

20  _ 21-Oct 

21  _ 21-Oct 

22  _ 21-Oct 

23  _ 23-Oct 

24  _ 23-Oct 

25  _ 24-Oct 

26  _ 25-Oct 

27  _ 25-Oct 

28  _ 25-Oct 

29  _ 25-Oct 

30  _ 31-Oct 

31  _ 31-Oct 

32  _ 1-Nov 

33  1-Nov 


Start  time 

(UTC) 

16:49 

19:47 

13:06 

16:42 

19:32 

12:22 


15:04 


22:25 


16:10 


15:44 


22:28 

21:19 


22:15 


End  time 
(UTC) 
17:20 
20:30 
13:32 
17:12 
19:53 
12:41 
15:31 
18:59 
13:23 
16:38 
19:04 
12:26 
15:42 
18:23 
22:34 
14:36 
17:27 
20:12 
11:21 
13:19 
14:43 
16:26 
13:18 
16:17 
12:28 
16:10 
18:51 
20:32 
22:45 
21:41 
22:33 
13:02 


Duration 

(Minutes) 

31 

j47 _ 

23 _ 

30 _ 

20 _ 

19 


#  of 

up/downs 

7 

1_ _ 

_5 _ 

1_ _ 

_6 _ 

7 


1-Nov 


Tables  3  and  4  catalog  all  MASL  measurements  from  the  RHIB/work  boat  during 
CASPER  pilot  and  CASPER  east.  The  tables  list  the  deployment  date,  start  and  end  time 
in  UTC,  duration  in  minutes  and  number  of  up/down  profiling  made  during  each  boat 
deployment.  Eigure  8  shows  an  example  of  the  MAPS  measured  profiles  for  02  May 
2015,  where  each  symbol  represents  the  measured  values  of  potential  temperature, 
specific  humidity,  and  the  derived  modified  refractivity  from  multiple  up/down  profiling 
in  one  location.  Although  scattered,  the  near  surface  gradients  in  potential  temperature 
and  M  in  the  lowest  few  meters  are  apparent  in  this  figure.  Extensive  data  processing  and 
mean  MASE  profile  retrieval  are  discussed  in  Chapter  IV. 


Theta  (K)  Sp.  hum  (g/kg)  M  (M  units) 

Eigure  8.  An  example  of  vertical  profiles  of  potential  temperature,  specific 
humidity  and  modified  refractivity  from  MAPS.  This  measurement  was 
made  on  02  May  2015  during  the  CASPER  pilot  experiment. 
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IV.  DATA  PROCESSING  AND  RESULTS 


A.  DATA  QUALITY  CONTROL 

We  performed  extensive  quality  control  on  the  MAPS-based  dataset.  This 
includes  identifying  errors  in  the  GPS  altitude,  deriving  pressure  altitude  as  altitude 
variables,  and  removal  of  erroneous  data  in  the  original  dataset.  This  section  will  outline 
these  efforts. 

1.  Altitude  Correction 

There  are  two  ways  of  obtaining  altitude  from  most  rawinsonde  sounding 
systems:  GPS  altitude  or  pressure  altitude.  GPS  altitude  has  a  known  global  average  error 
of  15  m  in  the  vertical  position  according  to  the  United  States  Department  of  Defense 
publication  on  GPS  performance  standard  (the  United  States  Department  of  Defense, 
2008),  which  is  due  to  the  design  of  the  GPS  system.  Expensive  GPS  receiver,  such  as 
the  RTK  or  Novatel  receivers  (https://en.wikipedia.org/wiki/Real_Time_Kinematic, 
http://www.novatel.com/an-introduction-to-gnss/chapter-5-resolving-errors/real-time- 
kinematic-rtk/)  have  an  accuracy  of  20  cm  or  less,  but  requires  a  base  station  similar  to 
differential  GPS.  Thus,  near  surface  high  accuracy  altitude  measurements  are  possible  but 
can  be  very  expensive.  The  high  end  GPS  receivers  are  not  available  in  the  iMet  sondes 
used  in  this  study 

Figure  9  shows  an  example,  taken  on  02  May  2015,  of  the  variation  of  GPS 
altitude  as  the  sensor  was  tethered  up  and  down.  It  is  seen  that  some  of  the  altitudes  were 
below  mean  sea  level,  which  was  inconsistent  with  visual  observations.  The 
corresponding  pressure  measurements  are  also  shown  in  Figure  9.  The  errors  in  GPS 
altitude  significantly  impact  the  vertical  profiles  of  scalars  as  seen  in  Figure  10.  Figure  10 
shows  temperature  and  relative  humidity  as  a  function  of  GPS  altitude  for  the  same 
measurements  in  Figure  9.  Here,  approximately  15-20%  of  the  data  points  were  shown 
below  mean  sea  level  as  indicated  by  the  shaded  area.  The  near  surface  profiles  also  show 
significant  scattering  at  these  ‘sub-surface’  levels. 
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21:30  21:32 


21:34  21:36 

Time  in  UTC 


Time  in  UTC 

Figure  9.  Pressure  and  GPS  altitude  from  02  May  2015  measurements. 


21:38  21:40 
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Figure  10.  Temperature  and  relative  humidity  as  a  function  of  GPS  altitude. 
The  measurements  were  made  on  02  May  2015.  The  shaded  area  denotes 
the  data  points  that  were  below  the  surface  by  GPS  altitude. 


To  mitigate  the  apparent  problems  with  GPS  altitude,  we  attempted  to  use  altitude 

derived  from  the  pressure  measurements,  the  so-called  ‘pressure  altitude’.  The  pressure 

altitude  was  derived  from  the  hypsometric  equation  using  pressure,  temperature,  and 

humidity  from  the  MAPS  system  as  input.  Surface  pressure  was  obtained  from  the  mast 

on  the  ship.  Figure  11  shows  the  same  variable  as  in  Figure  10,  except  with  the  derived 

pressure  altitude.  The  new  profiles  represent  significant  improvements  over  the  previous 
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one.  As  expected,  the  negative  altitudes  in  GPS  height  is  eliminated.  More  importantly, 
the  near  surface  data  points  show  significantly  less  scatter  compared  to  those  in  Figure 
10.  Apparently,  the  pressure  altitude  can  render  the  surface  layer  profile  characteristics 
much  better.  It  will  be  used  in  all  subsequent  analyses  as  default  altitude. 


Figure  11.  Temperature  and  relative  humidity  as  a  function  of  pressure  altitude. 

The  measurements  were  made  on  02  May  2015. 

2.  Erroneous  Data  Removal 

During  the  measurement  phase,  we  usually  turned  on  the  iMet  sonde  prior  to  the 
RHIB  leaving  the  dock  or  being  lowered  into  water  from  the  mother  ship.  Data  recording 
continued  throughout  the  cruise  until  the  RHIB  ship  came  back  to  dock  or  the  ship.  There 
are  inevitably  sections  of  data  outside  the  intended  measurements  period  such  as  transit 
time  or  sonde  up/down  measurement  setup  time.  This  erroneous  data  must  be  removed 
during  the  data  processing.  The  erroneous  data  may  occur  during  the  up/down  flight,  too, 
when  operators  adjusted  camera  or  the  tether  line.  Figure  12  shows  the  vertical  profiles  of 
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potential  temperature  and  specific  humidity  as  a  function  of  height  for  the  initially 
selected  period.  Figure  13  depicts  the  time  variation  of  the  pressure  altitude  showing  the 
up/down  sampling  of  the  tethered  balloon  based  sensor.  The  different  colors  in  both 
figures  represent  the  data  from  one  up/down  pair  during  this  set  of  measurements  at  the 
same  location.  Overall,  we  see  good  constancy  in  the  range  of  the  sampled  temperature 
and  humidity,  with  the  exception  of  the  blue  lines.  Figure  13  indicates  that  the  blue  line 
was  from  the  first  up/down  soundings  of  this  set.  It  is  very  likely  that  the  sensor  was  close 
to  the  operator  or  was  placed  somewhere  on  the  boat  initially  and  was  not  yet  in 
equilibrium  with  the  environment  during  the  first  sounding.  We  therefore  have  a  good 
reason  to  discard  this  portion  of  the  data.  The  green  line  also  shows  some  outliers.  Figure 
13  shows  the  sensor  moved  up  very  slowly  at  about  32  m  and  there  were  some  missing 
data  in  the  downward  sounding.  We  suspect  something  was  going  on  during  some  portion 
of  the  measurements.  The  grey  vertical  bars  in  Figure  13  indicate  the  periods  where  the 
data  was  removed. 
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Figure  12.  Potential  temperature  and  specific  humidity  as  a  function  of  height 
from  the  initially  selected  data  segment.  The  measurements  were  made  on 
31  October  2015  at  21 19  UTC  during  CASPER  East. 
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Figure  13.  Pressure  altitude  variation  of  the  initially  selected  data  segment  for 
the  same  balloon  flight  as  in  Figure  12.  The  shaded  area  denotes  the  data 
section  that  was  removed  in  the  final  quality  controlled  data  segment. 


Figure  14  shows  the  same  potential  temperature  and  specific  humidity  profiles 
after  the  erroneous  data  has  been  eliminated.  It  shows  a  tight  grouping  of  all  the  profiles 
and  shows  significantly  less  variation  than  in  Figure  12. 
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Figure  14.  Potential  temperature  and  specific  humidity  as  a  function  of  height 
after  erroneous  data  has  been  eliminated. 
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Figure  15  shows  another  example  taken  on  21  October  2015  at  1430  UTC  during 
CASPER  East  where  erroneous  data  are  present.  The  bad  data  occurred  mostly  at  the 
early  part  of  the  sampling,  especially  when  a  new  operator  was  onboard.  Figure  16  shows 
the  areas  in  grey  during  which  the  data  was  removed  from  the  final  dataset.  Finally, 
Figure  17  shows  a  cleaned  up  version  of  the  profile. 


Pot.  temperature  (K) 


Sp.  humidity  (g/kg) 


Figure  15.  Potential  temperature  and  specific  humidity  as  a  function  of  height 
from  the  originally  selected  data  section.  The  measurements  were  made  on 
21  October  2015  near  1430  UTC. 


Date  (UTC) 

Figure  16.  Pressure  altitude  variation  of  the  initially  selected  data  segment  for 
the  same  balloon  flight  as  in  Figure  15.  The  shaded  area  denotes  the  data 
section  that  was  removed  in  the  final  quality  controlled  data  segment. 
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Figure  17.  Potential  temperature  and  specific  humidity  as  a  function  of  height 
after  erroneous  data  has  been  eliminated. 


B.  GENERATING  MEAN  PROFILES 

The  MASL  is  inherently  turbulent.  Any  single  profiling  through  this  layer 
produce  a  single  measurement  at  each  level  and  hence  the  profile  from  such  measurement 
gives  a  composite  vertical  variation  seen  in  a  short  time  period  of  a  few  minutes  which  is 
the  time  to  complete  one  up  or  down  sounding.  With  multiple  up  and  down  soundings, 
the  variations  at  each  level  can  be  identified  as  seen  in  Figures  14  and  17.  This  section 
discusses  such  variability  and  the  methods  to  generate  a  mean  profile  from  each  set  of 
measurements  at  one  location. 

1.  Variability  in  the  Marine  Surface  Layer  Profile 

Figure  18  shows  the  potential  temperature,  specific  humidity  and  modified 
refractivity  as  a  function  of  height  for  a  single  profile.  The  observations  have  been  placed 
in  a  2  m  bin  averages  in  order  to  determine  the  mean  profile  for  all  three  variables 
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Figure  18.  Bin-averaged  profile  and  its  corresponding  variability  represented  by 
the  standard  deviation  of  the  data  point  in  each  vertical  bin.  Potential 
temperature,  specific  humidity,  and  the  modified  refractivity  are  shown  in 
this  figure.  The  measurements  were  made  on  15  October  at  1222  UTC. 


Figure  19  shows  bin-averaged  standard  deviations  of  potential  temperature  and 
specific  humidity  as  a  function  of  height  from  all  CASPER-East  MAPS  profiles.  The 
figure  shows  that  the  standard  deviation  above  10  m  is  in  general  less  than  0.3  °C  for 
potential  temperature  and  0.5  g  kg'^  for  specific  humidity,  although  the  standard 
deviations  varied  significantly  among  different  set  of  profiles.  However,  large  variations 
are  seen  in  the  lowest  10  m,  where  we  expect  large  vertical  gradient  based  on  MOST 
theory.  Overall,  the  standard  deviations  in  these  lowest  levels  are  less  than  0.5  “C  and  1.0 
g  kg'^  for  potential  temperature  and  specific  humidity,  respectively 
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Figure  19.  Bin-averaged  standard  deviations  of  potential  temperature  and 
specifie  humidity  as  a  function  of  height.  The  values  shown  are  composite 
from  all  CASPER  East  profiles. 


Eigure  20  separates  the  calculated  mean  standard  deviation  of  potential 
temperature  and  specific  humidity  for  all  vertical  bins  of  each  profiling  set  obtained 
during  CASPER  East.  These  vertically  averaged  variabilities  seem  to  be  rather  consistent 
from  day-to-day  and  at  different  locations,  except  for  the  last  few  profiles  obtained  over 
the  Gulf  Stream  region.  Eigure  21  shows  the  same  dataset  separated  by  date  of 
measurements.  One  can  see  that  the  variability  in  each  day  from  different  locations  has 
about  the  magnitude  as  the  day-to-day  variation. 
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Figure  20.  Mean  standard  deviation  of  potential  temperature  and  specifie 
humidity  from  each  MAPS  sampling  set  during  CASPER  East.  The 
horizontal  axis  is  the  sequential  sounding  number. 


Figure  21.  Same  as  in  Figure  20,  except  for  using  dates  in  local  time  as 

horizontal  axis. 


In  Figure  22,  we  examined  how  the  profile  mean  variability  differs  with  the 
number  of  up/down  profiles  (each  up/down  of  the  balloon  gives  two  profiles).  There 
appears  to  be  no  clear  trend  as  to  how  the  number  of  repeated  profiles  affect  the  sample 
variability.  In  other  words,  the  number  of  samples  we  used  for  these  measurements  seem 
to  be  a  sufficient  representation. 
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Figure  22.  Mean  standard  deviation  of  potential  temperature  and  specific 
humidity  from  each  MAPS  sampling  set  during  CASPER  East.  The 
horizontal  axis  is  the  number  of  profiles  taken  for  the  particular 

set  of  profiles. 


2.  Methods  to  Generate  Mean  Profile 

There  are  several  methods  developed  to  obtain  the  mean  profile  of  the  MASE 
from  each  MAPS  profiling  set.  The  first  method,  the  bin  averaged  method,  was  discussed 
in  the  previous  subsection.  We  used  2-m  vertical  bins  for  this  analysis.  The  size  of  the  bin 
gives  enough  samples  within  each  bin  and  was  small  enough  to  allow  the  distinguishable 
difference  in  adjacent  bins  near  the  surface  where  the  gradients  are  the  largest.  The 
second  method  uses  polynomial  fit  of  the  dataset  as  a  function  of  height.  Polynomial  fit  is 
a  curve  fitting  method  where  the  data  is  approximated  using  a  polynomial  function.  In 
this  study,  a  7th  order  polynomial  function  was  used  to  determine  the  mean  profiles.  This 
higher  order  fit  provides  more  flexibility  in  the  shape  of  the  profiles  and  allows  the 
representation  of  a  wide  range  of  vertical  variability.  The  third  method,  outlined  in  Kang 
and  Wang  (2016),  is  a  least-squares  optimization  method  that  utilizes  a  weighted  cost 
function  based  on  the  framework  of  MOST  with  the  cost  being  optimized  using  a  quasi- 
Newton  method.  The  method  assumes  the  samples  are  from  the  ensemble  of  profiles  that 
follow  MOST.  All  levels  of  the  profile  measurements  are  used  and  the  method  finds  the 
profile  that  minimizes  the  cost  function  with  a  specialized  weighting  function.  This 
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method  incorporates  all  of  the  measurements  in  the  derivations  of  the  profiles  and  fluxes. 
Details  of  this  method  are  documented  in  Kang  and  Wang  (2016). 

Figure  23  shows  the  mean  profiles  for  potential  temperature  and  specific  humidity 
taken  on  15  October  2015  generated  from  all  three  methods.  The  result  from  the 
polynomial  fit  (green  curve)  provides  a  smooth  curve  with  good  indication  of  the  mean 
gradients  in  the  low  levels.  For  the  most  part,  it  agrees  well  with  the  simple  bin  averaged 
profiles.  These  two  methods  are  based  on  the  measurements  only  with  no  assumptions. 
The  result  from  the  optimization  method  is  required  to  follow  the  Monin-Obukhov 
similarity  theory.  The  deviation  between  the  profile  from  the  optimization  method  and 
the  bin  averaged  or  poly  fit  methods  indicates  the  adequacy  of  MOST  theory  to  this 
profile.  The  results  in  Figure  23  show  that  most  of  the  differences  are  in  the  lowest  10  m 
of  the  surface  layer,  which  is  expected,  as  this  is  the  layer  of  the  largest  gradient.  One 
should  also  keep  in  mind  that  the  measurement  accuracy  is  likely  degraded  in  this  layer, 
because  of  the  strong  dependence  on  the  measured  quantities  of  height,  which  is  reflected 
in  the  larger  standard  deviations  in  the  lowest  levels  (Figures  18  and  19). 
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Figure  23.  An  example  of  mean  profiles  of  potential  temperature  and  specific 
humidity  generated  from  the  three  methods  discussed  in  the  text.  The  data 
was  obtained  on  15  October  2015  at  1222  UTC  (same  as  in  Figure  18). 
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Figure  24  shows  another  example  of  mean  generated  profiles  generated  from  all 
three  methods.  In  this  example,  all  three  methods  yielded  similar  potential  temperature 
profiles,  an  indication  that  the  Monin-Obukhov  Similarity  Theory  describes  this  case 
adequately.  As  in  the  example  in  Figure  23,  the  lowest  few  meters  does  not  show  as  much 
vertical  gradient  as  the  MOST  theory.  This  is  partly  due  to  the  limited  vertical  levels  in 
this  extremely  large  gradient  zone. 


Figure  24.  Mean  profiles  of  potential  temperature  and  specific  humidity  as  a 
function  of  height  taken  on  31  October  2015  at  21 19  UTC. 


C.  OBTAINING  TURBULENT  FLUXES 

The  golden  standard  method  to  obtain  fluxes  is  the  eddy  correlation  method, 
which  requires  high-rate  measurements  of  the  temperature,  humidity,  and  wind 
components.  This  type  of  data  was  available  from  the  R/V  Hugh  Sharp  and  the  R/V 
Atlantic  Explorer  during  CASPER  East.  The  MAPS  measurements  from  the  R/V  Sharp’s 
work  boat  was  always  very  close  to  the  R/V  Sharp,  we  therefore  use  the  derived  flux 
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from  the  high-rate  measurements  on  the  R/V  Sharp  to  generate  the  true  surface  fluxes  of 
heat,  water  vapor,  and  momentum. 

Turbulent  fluxes  can  also  be  calculated  using  surface  flux  parameterization  such  as 
the  COARE  algorithm.  The  input  mean  quantities  can  be  obtained  from  the  ship  mast  on 
the  RfV  Sharp.  We  can  also  use  the  scalar  quantities  measured  by  the  MAPS  as  input  to  the 
COARE  algorithm  to  obtain  the  parameterized  surface  fluxes.  A  comparison  of  the  fluxes 
with  all  three  methods  are  shown  in  Eigure  25.  To  calculate  the  fluxes  from  the  mast  we 
used  temperature,  humidity  and  wind  values  obtained  from  the  ships  mast  and  input  these 
values  into  the  COARE  algorithm  to  obtain  the  parameterized  flux  values.  To  obtain  the 
fitted  profile  fluxes  we  used  the  mean  values  from  our  profile  observations  with  the 
exception  of  winds.  Since  the  sonde  is  tethered,  no  measurements  of  winds  were  available 
from  the  sonde  measurements.  As  a  result,  the  12  m  winds  from  the  ships  mast  were  used. 
Hence,  the  parameterized  fluxes  from  the  mast  and  the  MAPS  use  the  same  wind  input, 
which  explains  some  of  the  similarities  for  fluxes  from  these  methods  (Eigure  25). 


40 


Oct  13  Oct  16  Oct  19  Oct  22  Oct  25  Oct  28  Oct  31 

Date  (LT) 

Figure  25.  Momentum  (MF),  latent  heat  (LHF)  and  sensible  heat  flux  (SHF) 
during  each  CASPER  East  profiling  period. 
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D.  EVAPORATIVE  DUCT  HEIGHT  AND  STRENGTH 


During  CASPER  East,  there  were  36  sets  of  profiling  data  collected  at  different 
times  and  locations.  We  are  able  to  obtain  the  mean  profiles  from  each  measurement  set 
and  obtain  the  evaporation  duct  information  from  the  M-profiles.  Eive  examples  of  the 
derived  evaporation  duct  and  the  corresponding  profiles  are  shown  in  Eigures  26-30. 
These  cases  were  chosen  based  on  the  surface  layer  thermal  stability  represented  by  air- 
sea  temperature  differences,  which  represent  strongly  stable,  stable,  nearly  neutral, 
unstable,  and  strongly  unstable  thermal  stability.  The  blue  dots  are  the  observed  data 
points,  the  red  lines  are  the  polynomial  fitted  mean  profiles,  and  the  green  asterisk  (*)  on 
the  modified  refractivity  plot  represents  the  evaporative  duct  height.  The  value  of  EDH 
from  each  case  is  given  at  the  top  of  the  figures  together  with  the  time  in 
yyyymmddhhmm  (UTC)  and  the  air-sea  temperature  difference  for  this  case.  The  SST 
was  obtained  by  the  high  accuracy  integrated  infrared  SST  autonomous  radiometer 
(ISAR)  on  the  R/V  Sharp.  The  air  temperature  was  obtained  from  the  ship  mast.  Since  the 
workboat  was  always  in  the  vicinity  of  the  R/V  Sharp,  the  air-sea  temperature  difference 
should  be  representative  of  the  work  boat  environment  as  well. 


strongly  stable 

AT-SST=1.6°C  201510131947  EDH  =43.8  m 


295  296  297  298  10  12  14  16  345  350  355  360 


Pot.  temperature  (K)  Sp.  hum  (g/kg)  Mod.  refractivity  (M  units) 

Eigure  26.  Mean  profile  of  potential  temperature,  specific  humidity  and 
modified  refractivity  as  a  function  of  height  for  strongly  stable  case  obtained 
on  13  October  2015  at  1947  UTC. 
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Figure  27.  Mean  profile  of  potential  temperature,  speeifie  humidity  and 
modified  refraetivity  as  a  funetion  of  height  for  a  ease  of  stable  stratification 
measured  on  25  October  2015  at  1811  UTC. 
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Figure  28.  Mean  profile  of  potential  temperature,  specific  humidity  and 
modified  refractivity  as  a  function  of  height  for  a  near  neutral  case. 
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Unstable 

AT-SST= -1.4°C  201510201657 


Pot.  temperature  (K)  Sp.  hum  (g/kg)  Mod.  refractivity  (M  units) 

Figure  29.  Mean  profile  of  potential  temperature,  specific  humidity  and 
modified  refractivity  as  a  function  of  height  for  an  unstable  case. 
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Figure  30.  Mean  profile  of  potential  temperature,  specific  humidity  and 
modified  refractivity  as  a  function  of  height  for  a  strongly  unstable  case. 


Figures  26-30  show  a  substantial  range  of  EDH  from  the  five  cases.  The  deepest 
EDH  was  found  in  the  strongly  stable  case  followed  by  the  moderately  stable  case.  In  the 
strongly  unstable  case,  no  evaporative  duct  height  was  identified,  this  is  possibly  due  to 
the  strong  mixing  in  the  surface  layer  caused  by  large  buoyant  eddies. 

Eigures  31  shows  a  comparison  of  EDH  and  EDS  as  a  function  of  time.  The 
observed  EDH  and  EDS  are  directly  from  the  mean  profiles  obtained  using  the 
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polynomial  fit  method.  The  modeled  EDH  and  EDS  use  three  set  of  input  variables  from 
the  bow  mast  of  the  R/V  Sharp,  the  small  mast  on  the  workboat,  and  from  the  mean 
MAPS  profile  at  a  given  (12  m)  height.  The  data  in  generally  is  in  good  agreement  but  we 
do  see  some  inconsistencies  within  the  data.  Apparently,  the  modeled  EDS  using  the  ship 
bow  mast  data  as  an  input  consistently  overestimates  EDS. 


Eigure  31.  Comparison  of  the  measured  and  model  evaporation  duct  properties 
for  the  entire  CASPER-East  profile  set. 


50 


40 


-^30 

I 

Q 

la  20 


10 


O  Poly  fit 

O  Ship  mast-12  m  data 
O  Boat  mast-1  m  data 
O  Poly  fit-12m  data 


-4  -2  0 

Tair-SST  (“C) 


80 


70 


_^60 

t/T 

1  50 


W  40 

Q 

LU 

30 


20 


10 

-6  -4  -2  0  2  4 

Tair-SST  (“C) 


< 

0 

o 

n  < 

°o 

!pcP^ 

< 

o 

o  < 

A 

Q  0 

% 

O  /■ 

o 

o 

Ob 

8 

0  ^ 

0 

“a 

Eigure  32.  Observed  or  modeled  evaporation  duct  properties  as  a  function  of 
ASTD.  Input  data  to  the  model  was  obtained  from  R/V  Sharp’s  bow  mast  at 
12  m,  1  m  boat  mast,  and  12  m  polynomial  fit  over  static  stability  criteria  for 
the  entire  CASPER-East  profile  set. 


45 


Figure  32  shows  the  same  evaporation  duct  properties  as  in  Figure  31,  except  as  a 
function  of  air-sea  temperature  difference  (ASTD).  Negative  ASTD  denotes  unstable 
stratification,  positive  for  stable  stratification.  Here,  we  observe  small  variability  of  the 
EDH  in  the  unstable  conditions  and  strong  sensitivity  in  stable  conditions.  This  result  is 
consistent  with  the  findings  of  Cherrett  (2015)  using  the  COARE  algorithm.  However, 
there  was  still  a  lack  of  observational  data  in  the  stable  regime.  The  measured  EDS 
follows  the  MOST  calculated  values  well  and  indicated  weak  ducts  in  stable  conditions 
and  large  enhanced  EDS  in  the  unstable  conditions.  In  Eigure  33  where  the  comparisons 
are  made  for  cases  from  the  west  coast  during  CASPER  Pilot,  the  model  seems  to 
overestimate  EDH. 


Eigure  33.  Measured  evaporative  duct  height  verses  estimated  using  COARE 

algorithm. 
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V.  SUMMARY  AND  CONCLUSIONS 


A.  NPS  MAPS  CONCLUSIONS 

Accurate  depiction  of  the  temperature,  humidity,  and  wind  profiles  in  the  MASL 
are  critical  in  quantifying  the  effects  of  the  atmosphere  on  electromagnetic  wave 
propagation  and  in  particular  evaporative  ducts.  However,  the  vertical  profiles  in  the 
lowest  50  m  of  the  atmosphere  have  not  been  well  measured  because  of  the  impact  of  the 
large  ships  or  platforms  on  their  immediate  environment.  NPS  has  developed  a  Marine 
Atmospheric  Profiling  System  (MAPS)  that  is  capable  of  making  repeated  measurements 
of  the  MASL  in  a  relatively  undisturbed  environment.  The  MAPS  is  comprised  of  a  MET 
mast,  tethered  balloon,  a  radiosonde  and  receiving  and  data  display  package  along  with  a 
sea  surface  thermistor.  The  system  is  designed  to  obtain  surface  layer  profiles  as  well  as 
one  level  of  wind  for  data  retrieval  and  analysis  later  on.  The  system  prototype  was  tested 
in  CAPER-Pilot  and  improved  in  CASPER-East,  15  data  sets  were  made  during 
CASPER-Pilot  and  36  during  CASPER-East. 

Careful  data  quality  control  was  applied  to  the  MAPS  profiles.  Issues  with  GPS 
altitude  were  identified  that  may  introduce  the  structure  and  characteristics  of  the  profiles 
considerably.  This  problem  was  resolved  by  using  pressure  altitude.  Erroneous  data, 
mostly  in  the  beginning  of  the  data  collection  at  each  location,  was  removed  by  using 
prior  knowledge  of  the  problem  domain  and  by  using  common  sense.  The  variability  of 
the  MASE  temperature  and  humidity  at  each  level  was  clearly  seen  from  multiple 
up/down  measurements  at  a  single  location.  These  variations  are  result  of  the  natural 
variability  in  a  turbulent  environment,  which  points  out  to  potentially  large  uncertainties 
in  using  a  single  vertical  profile  to  represent  the  mean  atmospheric  conditions.  The  mean 
profile  and  its  variability  was  obtained  by  sectioning  the  vertical  dimension  into  2-m 
depth  bins  and  calculate  the  mean  and  variance  in  each  bin.  The  bin-averaged  vertical 
profiles  clearly  showed  the  near  surface  gradients  in  both  temperature  and  humidity, 
which  allowed  us  to  identify  the  evaporation  duct  properties  from  the  derived  mean 
modified  refractivity  profile.  The  variances  of  potential  temperature  and  specific 
humidity  did  not  change  significantly  with  height  above  the  surface  and  are  consistent 
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for  most  of  the  profiles  with  the  exception  of  profiles  taken  over  the  Gulf  Stream  on  01 
November  2015.  The  variances  also  did  not  vary  significantly  for  different  number  of 
sample  profiles.  In  addition  to  the  simple  bin  averaging  approach,  a  polynomial  fit  of  the 
data  points  were  also  attempted  to  generate  the  ‘mean’  profile.  The  7*  order  polynomial 
fit  and  the  bin  average  method  gave  rather  consistent  profiles.  In  addition,  Kang  and 
Wang  (2016)  introduced  a  least- square  optimization  method  that  fit  the  data  points  to  a 
profile  that  followed  the  MOST  in  the  surface  layer.  This  method  was  also  applied  to  the 
CASPER  East  and  some  of  the  CASPER  Pilot  profiles.  Good  comparison  between  this 
method  and  the  purely  observation  based  profiles  indicate  the  validity  of  the  MOST 
theory,  which  is  the  case  in  some  of  the  profiles.  However,  there  are  also  profiles  that 
showed  significant  deviations. 

Using  data  from  a  specific  level  on  the  MAPS  profiles  as  an  input  to  surface  flux 
parameterization  algorithms,  one  can  obtain  surface  fluxes  of  sensible  and  latent  heat. 
Eluxes  generated  this  way  compared  well  with  those  from  the  eddy  correlation  method  or 
the  parameterized  fluxes  using  measurements  from  the  bow  mast  on  the  R/V  Sharp  near 
by  the  workboat  during  CASPER  East.  Since  the  wind  speed  used  in  the  parameterized 
flux  calculation  was  the  same,  the  good  comparison  of  the  parameterization  essentially 
reveals  that  the  mean  temperature  and  specific  humidity  from  the  ship  bow  mast  and  the 
MAPS  profile  at  the  bow  mast  level  are  consistent.  The  good  comparison  between  the 
eddy  covariance  and  the  parameterized  fluxes  suggest  the  validity  of  the  COARE  flux 
algorithm,  which  was  used  in  calculating  fluxes  from  the  mean  measurements. 

The  main  objective  of  developing  the  MAPS  is  to  obtain  directly  measured 
evaporation  duct  properties.  Using  the  surface  layer  profiles  from  CASPER  East  and 
CASPER  Pilot  field  campaigns,  we  were  able  to  obtain  evaporation  ducts  properties  in 
different  stability  and  wind  conditions  shown  in  five  examples  in  this  thesis.  Shallow  but 
strong  ducts  were  found  in  the  unstable  surface  layers  and  deep  but  weak  ducts  were 
found  in  the  stable  surface  layers.  The  measured  evaporation  duct  properties  can  be  used 
to  evaluate  evaporation  duct  models.  In  this  thesis,  we  used  the  modified  COARE 
algorithm  to  derived  the  modeled  evaporation  ducts.  A  Comparison  of  measured 
evaporative  duct  height  with  COARE  surface  layer  model  estimated  evaporative  duct 
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height  from  CASPER  pilot  experiment  showed  that  model  generally  overestimated  the 
duct  heights  in  cases  during  CASPER  pilot.  No  persistent  trend  of  the  model-observation 
comparison  was  observed  for  cases  obtained  during  CASPER  East.  Eurther  in-depth 
analyses  are  needed  with  this  dataset  to  identify  the  conditions  where  evaporation  models 
need  further  improvements. 

B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  NPS  MAPS  is  fully  developed  and  functions  very  well  in  CASPER  pilot 
testing  and  during  CASPER  East.  We  have  also  developed  methods  to  obtain  mean 
temperature  and  humidity  profiles  for  calculation  of  the  mean  evaporation  duct  height. 
Euture  work  should  be  associated  with  the  use  of  the  system  and  deploying  it  in  different 
conditions  so  that  more  data  can  contribute  to  a  better  understanding  of  the  evaporative 
duct  modeling. 
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